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ABSTRACT 

We present results of a measurement of the shape of the density profile of galaxy groups 
and clusters traced by 43 335 Luminous Red Galaxies (LRGs) with spectroscopic 
redshifts from the Sloan Digital Sky Survey (SDSS). The galaxies are selected such that 
they are the brightest within a cylindrical aperture, split into two luminosity samples, 
and modeled as the sum of stellar and dark matter components. We present a detailed 
investigation of many possible systematic effects that could contaminate our signal and 
develop methods to remove them, including a detected intrinsic alignment for galaxies 
within 100/i~^kpc of LRGs which we remove using photometric redshift information. 
The resulting lensing signal is consistent with NFW profile dark matter halos; the SIS 
profile is ruled out at the 96 (conservatively) and 99.96 per cent confidence level (CL) 
for the fainter and brighter lens samples (respectively) when we fit using lensing data 
between 40 /i~^kpc and 2 h~^Mpc with total signal-to-noise of 19 and 25 for the two 
lens samples. The lensing signal amplitude suggests that the faint and bright sample 
galaxies typically reside in haloes of mass (2.9 ± 0.4) x 10^^h~^MQ and (6.7 ± 0.8) x 
10^^h~^MQ respectively, in good agreement with predictions based on halo spatial 
density with normalization lower than the 'concordance' tJg = 0.9. When fitting for 
the concentration parameter in the NFW profile, we find c — 5.0 ± 0.6(stat) ± l(sys), 
and c = 5.6 ± 0.6(stat) ± l(sys) for the faint and bright samples, consistent with 
ACDM simulations. We also split the bright sample further to determine masses and 
concentrations for cluster-mass halos, finding mass (1.3 ± 0.2) x lO^'^h'^ Mq for the 
sample of LRGs brighter than —22.6 in r. We establish that on average there is a 
correlation between the halo mass and central galaxy luminosity relation that scales 
as M oc L'^. 

Key words: gravitational lensing — galaxies: haloes. 



1 INTRODUCTION 

Elliptical galaxies, particularly Luminous Red Galaxies 
(LRGs), are good tracers of the most massive galactic ha- 
los. They have a number of useful properties, including sim- 
ple, well-understood spectral energy distributions that can 
be modeled as a passively-evolving burst of star formation, 
with a constant low level of ongoing star for mation; variation 
with environment t hat is quantifi ed as well ( Eiscnst ein et all 
I2OOII: iBruzual fc Charlot..2003, : .Eisenstein et al...200ll : uni- 
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form photometric properties, including the fundamen- 
tal plane relat ing luminos i ty, radius, and central veloc- 
ity dispersion l|Fab ei]|l973t [Visvanathan fc Sandagd[l977| : 



DiorgQvski fc DavisI Il987l: iKormendy fc Diorgovskil 198! 
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Bower et alJll992l : iRoberts fc Havnesl ll994': Bcrnardi et al 
12003) of which the Faber- Jackson relation ( Fabcr & Jackso, 
Il97(fl is a projection; and simple light profiles, the de Vau- 
couleurs profile, with I oc exp(— r^''*). Galaxies on the mas- 
sive end of the red sequen ce contain the majority of the stel- 
lar mass of the universe (Fukugita ct al. 1998; 'Hogg et ajj 
^002; Bell ct al. 2004; Baldrv ct al. 2004; Cool ct al. 200^^ 
and hence their properties are of great interest. It was shown 
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that while eUiptical galaxies dominate over spirals in over- 
dense environments, these galaxies do reside in a wide va- 
riety of environments, f r om the field to the richest clus- 
ters (Morgan ct al." 1975; "A lbert et aPllQTTt iPonman et alJ 
\l99A: Vikhlinin ct al. .1999 ,: Loh & Strauss, 2006, in prep.), 
with the majority of the LRGs located in group to small 
cluster sized halos. There have also been ex tensive studies 
of LRG clustering Eisenstein et al.ll2005al lbt iMasiedi et alJ 



20051: IZehavi et alJl2005^ from 10 fe~^kpc to 100 h~^Mp 
scales. 

In addition to the extensive photometric studies of 
elliptical galaxies already in the literature, we would like 
to add a good understanding of their density profiles. 
This information is useful for several purposes. First, 
the dark matter profile can give us information that can 
be related to the formation history and merger history 
of these galaxies. Second, the shape of dark matter 
(DM) profiles can be pred i cted using N-body si mula- 
tions jNavarro et"ai] Il996t iFukushiee fc Makind Il997l: 



Kravtsov et alJ Il997l:' Moore et alJ 11993: lAvila-Reese et al 



19991 : iGhigna et al.ll2000t Ijing fc Sutdl2000l: iBuUock et al 

2001; 'Klvpin et al.' '200l'; 'Pukushigc fc Makin g l200lt 
Wechsler et al 2002: .Fukushige fc Making 2003: Zhao et all 



20031: iTasitsiomi et alJl2004l : Diemand et al.ll2005^ ■ though 
with some disagreement in the value of inner asymptotic 
logarithmic slope, possibly due to resolution issues, and 
hence LRGs in massive halos provide a particularly good 
way of testing the density profile predictions of ACDM. The 
brightest red galaxies have been shown to have small baryon 
fractions compared to fainter red galaxies and compared 
to typical blue galaxies ('Mandclbaum ct al. 2005d), so the 
effects of baryons on the DM profile are relatively small in 
this class of galaxies. Furthermore, since the virial radii are 
so large, and the concentration parameters decrease with 
mass, any changes in logarithmic slope around the scale 
radius occur at large scales relative to smaller galaxies, and 
hence are easier to measure using the tools available. 

While the matter distributions in elliptical galax- 
ies can be simply probed on small sc ales (within a few 
re) u s ing central velocity d i spersions jKronawitter e t alJ 
I2OOOI: Ide Zeeuw et al.l 120021: IPadmanabhan et all |2004) . 
finding probes of the dark matter profiles on larger scales 
is more challenging. Kinematic tracers such as satellite 
galaxies can give information out to tens of /i~^kpc. 
Hydrostatic analyses of X-ray intensity profiles, and strong- 
and weak-lensing constraints for individual clusters, are 
numerous but thus far do not g jive a single clear picture 
jAUed Il9 98>: "Clowe et al.' "2000"; Clowe 
Rusin fc Ma 2001; Sheldon ct al 



Schneideij 

2001; Athreva et alJ 



r picture 
i eil l200lt 



^^^^^^ , - '^Y^ - - 

20021: lAr abadiis et a l.i l2002t iKoopmans fc Treul l2003t 



Huterer fc Mai i2004»: IFerreras et al. 
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20051: I ^ [200^ iHumpln'e^^^jl I2OO6- 
Koopni£m^^^Lri200^ . Weak lensing is another tool 



that can be used out to scales of several /i~^Mpc, and 
the SDSS provides a particularly advantageous dataset 
due to its significant sky coverage and spectroscopic 
redshifts (which allow us to deter mine the profil e as a 
function of trans verse separation). iHoekstra et al.l ll200lll 
and Parker et al] |2005) have done a group lensing analysis 
with a much smaller sample (~ 100) groups from CNOC2; 
however, these small samples are not sufficient to allow a 
detailed determination of the shape of the density profile. 



Measurements of the 3-dimensional density profiles can 
be used to learn about cosmology as well. For example, the 
concentration parameter of the profile can be related to the 
halo mass via the dark matter power spectrum normaliza- 
tion as and the matter density f2,„. Thus, measuring the 
profile for several samples of varying halo mass can teach us 
about the underlying cosmology. 

We begin in S|5|with a description of the theory behind 
the measurement we are making. [0 has a full description 
of the dataset used, including the properties of the lens and 
source samples, and possible causes of systematic error. We 
present results of profile measurements in [21 ^^nd conclude 
in !gl 

Here we note the cosmological model and units used in 
this paper. All computations assume a flat ACDM universe 
with Qrn = 0.3, Qa = 0.7, and as ~ 0.9. Distances quoted 
for transverse lens-source separation are comoving (rather 
than physical) /i~^kpc, where Hg — 100 h kms~^Mpc~^. 
Likewise, AE is computed using the expression for E^^^ in 
comoving coordinates, Eq.|3] In the units used, Hq scales out 
of everything, so our results are independent of this quan- 
tity. Finally, 2-dimensional separations are indicated with 
the capital R, 3-dimensional radii with lower-case r (occa- 
sionally r may denote r-band magnitude as well). 



2 THEORY 

Here we describe the theory behind our attempts to mea- 
sure the three-dimensional density profile using the galaxy- 
galaxy weak lensing signal. 



2.1 Lensing 

Galaxy-galaxy weak lensing provides a simple way to probe 
the connection between galaxies and matter via their cross- 
correlation function 



(1) 



where 5g and 5m are overdensities of galaxies and matter, 
respectively. This cross-correlation can be related to the pro- 
jected surface density 



E(J?) 



(2) 



(where — R^ + x'^) which is then related to the observable 
quantity for lensing. 



AE(i?) = 7t(i?)E, = E(< 7?) - E(7?), 



(3) 



where the second relation is true only in the weak lensing 
limit, for a matter distribution that is axisymmetric along 
the line of sight (which is naturally achieved by our pro- 
cedure of stacking thousands of lenses to determine their 
average lensing signal). This observable quantity can be ex- 
pressed as the product of two factors, a tangential shear 7t 
and a geometric factor 



Ds 



(4) 



where Dl and Ds are angular diameter distances to the lens 
and source, Dls is the angular diameter distance between 
the lens and source, and the factor of (1 + zz,)~^ arises due 
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to our use of comoving coordinates. For a given lens redshift, 
E"'^ rises from zero at Zs = zl to an asymptotic value at 
Zs 3> Zl; that asymptotic value is an increasing function of 
lens redshift. 

For this paper, we are primarily interested in the contri- 
bution to the galaxy-mass cross-correlation from the galaxy 
halo profile itself (central galaxy Poisson term) , rather than 
from neighboring halos (halo-halo term), and hence 



E(i?) 



(5) 



The halo-halo term for galaxies in host halos can be modeled 
simply usi ng the galaxy-dark matter c ross-power spectrum 
as in, e.g., iMandelbaum et al ] i2005b), and is only impor- 
tant for R > 2 h~^Mpc, which is the maximum scale probed 
in this paper. 

2.2 Profiles 



In previous analyses (IMandelbaum et alJl2005blld|) . we have 
modeled the lensing signal as the sum of contributions from 
central galaxies (those that lie in host halos) and from satel- 
lite galaxies (those that lie in subhalos). For this work, to 
simplify interpretation by avoiding issues such as tidal strip- 
ping of satellites and the radial distribution of satellites 
within host halos, we use methods (described in H3.1II that 
enable us to isolate a sample of central galaxies with negli- 
gible contamination from satellites. Thus, for this work, our 
interpretation is much simpler, and can be related directly 
to p{r) for the lens sample under consideration. 

While more complicated models may be useful for tech- 
niques capable of resolving smaller scale information than 
this work (minimum transverse separation of 20 /i~^kpc), 
we model each galaxy as the sum of two components. The 
first component is a dark matter profile, with the following 
density profile: 



p{r) — _^ 

(r/rs) " (1 + r/r, 

which gives a logarithmic slope of 



d In p 
d In r 



= a-{a + f3) 



r/rs 



1 + r/rs 



(6) 



(7) 



or a for r <ti rs and — /3 for r . This profile wit h a = — 1 
and /3 = 3 is the NFW profile ijNavarro et al.ll99fi) . and this 
shape (with varying values of a) is a generic prediction of 
ACDM N-body simulations over a large range of masses. In 
addition to the logarithmic slopes, it can be defined by two 
parameters, c = r^ir/rs and M. The virial radius r^ir and 
ps can be related to M via consistency relations. The first 
is that the virial radius is that within which the average 
density is equal to 180p: 



47r 
M = — r. 
3 



(180p) 



(8) 



Note that this definition differs from the oft-used 200pcrit 
definition of mass by roughly 30 per cent for typical values 
of concentration. The second relation, used to determine ps 
from M and c, is simply that the volume integral of the den- 
sity profile out to the virial radius must equal the virial mass 
(though when computing the signal, we do not truncate the 
profiles beyond r^ir)- In principle, the NFW c is a weakly de- 
creasing function of halo mass, with a variety of expressions 



used te uUock et al ]|200ll : lEke et alJl2001^ , making this pro- 
file a one-parameter family of profiles, but in order to give 
generality to the fits, we allow it to be a free parameter. 

To generate the profile at a particular redshift, we start 
with a given mass, use Eq.|H|to get the comoving virial radius 
using the comoving matter density, so that M, a, /3, and c 
define the profile uniquely. When computing the lensing sig- 
nal, we are actually averaging over a range of redshifts, which 
means that since we define the profile in comoving coordi- 
nates, we expect that the measured profile will be somewhat 
blurred due to evolution of concentration with redshift. 

While the more general ellipsoidal versions of these den- 
sity profiles are seen universally in ACDM simulations, here 
we work only with the spherical versions, because our pro- 
cess of averaging over thousands of lens galaxies makes the 
results insensitive to halo triaxiality. This statement is not 
made on the basis of analytical proofs (and surely must be 
wrong if one requires an extremely high level of precision), 
but rather on the basis of the lensing signal measured in 
N-body simulations as an average over randomly-o riented 
triaxial NFW profiles. In lMandelbaum et alj (l2005d) . it was 
shown that spherical NFW profiles do an excellent job at 
describing the lensing signal for a variety of luminosity bins, 
with the masses and concentrations of the best-fit profiles 
related to the real masses and concentrations in the simu- 
lations in a particular way (to be discussed further below). 
The fit values were very good when the errorbars used on 
the signal were one-tenth of our current errorbars, so hence- 
forth we consider only spherical NFW fits rather than trying 
to account for the averaging of triaxial halos. 

The second component to the full p(r) that we con- 
sider is the baryon density. This includes not only the stel- 
lar mass, but also any baryons that did not form stars. 
iRead fc Trenthaml (120051) estimate that 80 per cent of the 
baryonic mass in galaxies is in the form of stars, though for 
elliptical galaxies this number is closer to 100 per cent, so 
there is very little gas that has not been converted to stars. 
As is commonly done for ellipticals, we model the p(r) for 
the stellar component as a Hernquist profile since it is a de- 
projected de Vaucouleurs profile. The Hernquist profile can 
be expressed via Eq.|H]with a = —1 and /3 = 4 (and hence 
falls off faster than an NFW profile on large scales). In the 
absence of further information about any "dark" baryonic 
component, we do not model its profile separately, with the 
expectation that the need for a dark baryonic component 
will be shown by poor fits in the inner regions. 

Finally, we cannot in reality assume that these density 
profiles are independent. The gravitational potential due to 
one will necessarily affect the form of the other, so the dark 
matter profile will not be expected to conform exactly to 
the ones seen in N-body simulations on small scales where 
baryonic matter is sig nificant. This effect is known as adia- 
batic contraction (ACjBlumentha^^al .'l986^.'G nedin et alJ 
^004, Gncdin 2005!, ISellwood fc McGa ugh 200|) and has 
been studied with analytic models and simulations. To 
model a d iabat ic contraction, we use the CONTRA^ code from 
iGnedinI (|200^ to obtain adiabatically-contracted density 
profiles p(r). Generally, on the scales we consider, which are 
several times the de Vaucouleurs radius, we will show that 

^ \http://\vww-astionomy.mps.ohio-state.edu/~ognedm/contra/\ 
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Figure 1. Density profile p{r) for tlie typical LRG described 
in the text, including the profile separately for dark matter and 
for baryons, the total profile, and the total profile including AC 
(barely distinguishable from the profile without AC), as labelled 
on the plot. 



Figure 2. Surface density and projected surface density 

AS(i?) for the typical LRG described in the text, including the 
profile separately for dark matter and for baryons, and the total 
profile without AC (AC produces a negligible difference at the 
relevant transverse separation). A vertical line at 40 /i~^kpcshows 
the minimum scale used for the fits. 



AC leads to an immeasurably small change in the profile, 
and lience is not necessary for this work. This is because of 
low baryonic fractions in these galaxies. 

Figure Q shows the density profile p(r) for a typical 
expected LRG profile, with Mnfw = 6 x IQ^^H'^Mq, 
Ms = 8 X lO"/i"^M0, c = 7, baryonic = 10 /i^^kpc, 
Q = -1, /3 = 3, E;:^ = 2/i"^ X 10"''pc^Mq\ with and with- 
out AC. The source for these numbers are the predicted LRG 
masses from the halo mass function, a calculation that will 
be described shortly. 

Given that we measure the lensing signal for projected 
separations i? > 20 /i~^kpc, and fit using i? > 40 /i~^kpc, 
one might wonder why a baryonic component is necessary 
for our modeling at all, given the values of the density pro- 
files in Fig. at these separations. We present Fig. |5| to 
answer this question. The top panel of this figure shows the 
projected surface density E as a function of transverse sep- 
aration for the same density profile as in Fig. and the 
bottom panel shows the surface density contrast, AE, which 
is the quantity of interest for lensing. As shown in the top 
panel of Fig. |5| the baryonic component produces a signif- 
icant increase in the surface density for transverse separa- 
tions R <~ 20 /i~^kpc. However, what is important for the 
computation of AE is that also leads to a significant steep- 
emng of E on these scales. This steepening of E leads to a 
noticable increase in AE even on larger scales due to the de- 
pendence on AE(< R), out to 100 /i~^kpc for this 'typical' 
profile. This effect is the reason why our modeling requires 
a baryonic component in AE for 20 < R < 100 /i~^kpc de- 
spite its relatively small contribution to p for similar scales. 

As mentioned above, the full lensing signal for host 
galaxies (Eq. |5J is the sum of the terms due to the cross- 
correlation between a host galaxy and its own halo (Eq. (SJ, 



and the halo-halo term (cross-correlation with the dark mat- 
ter in other halos). We must still cons ider the halo- halo term. 
Simple modeling of the h-h term as in ISeliakI J2000l) indicates 
that below 2 h~^Mpc, this term is <C AE due to the central 
profile. 

One factor that we must consider when relating the 
best-fit profile parameters to reality is the effect of aver- 
aging over a distribution of halo masses. We start with the 
halo mass function, dn/dM (which depends on redshift), to 
specify the distribution of halo masses: 

dn 



-dM 



-f{v)dv 



(9) 



dM M- 

where p is the mean matter density of the universe. The 
function f{i/) can we written in a universal form when de- 
fined in terms of peak height 



(10) 



a{M) 



where 5c = 1.686 is the linear overdensity at which a spher- 
ical perturbation collapses, and a{M) is the rms fluctua- 
tion in spheres containing on average mass M at an ini- 
tial time, extrapolated using linear theory to z. We use the 
ISheth fc TormenI (|W9^ version of the mass function. 



(11) 



wh ere v' = av, with a = .707 and p = 0.3. The origi- 
nal ^^s^^^checMe^ il974l) mass function corresponds to 
a = 1 and p = 0. The constant A is determined by mass 
conservation, requiring 



f(v)dV : 



(12) 







Once we have normalized the mass function, we can 
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compare the comoving density of the LRG sample to the 
predicted density of halos with mass above some minimum 
mass Mmin to estimate minimum halo masses, assuming 
that the majority of these massive halos are populated by 
red rather than blue galaxies. For a comoving density of 
riLRG = 1-2 X 10"''(/i/Mpc)^, as for the fuU LRG sample, 
the characteristic sample redshift suggests that 



-jTjdM = ulrg 
aM 



(13) 



for Mmin ~ lO^^ "* = 2.5 X 10"/i~^Mq. As will be shown, we 
also split the sample into two, such that the fainter sample 
contains roughly 2/3 of the galaxies. We can model this split 
by saying the fainter luminosity sample has the same Mmin 
and includes all halos up to some Mi (i.e., a top-hat filter 
on n(M)), using the fact that light roughly traces mass for 
these galaxies. Then, we find Mi using 

dn 2 

-dM = -ULRG (14) 



dM 



which gives A/i 



10^ 



= 6 X lO"/i~^M0. This Mi is 



then the minimum halo mass for the brighter LRG sample 
containing 1/3 of the LRGs. 

For these two samples containing a range of halo masses, 
we must then consider the relationship between AE aver- 
aged over n{M), i.e. 



(AE) 



/AE^dM 



(15) 



over the relevant mass range {Mmin to Mi, or Mi to 
oo). Because of difi'erences in expected concentration param- 
eters, one might expect that this averaging process will 
change not only the signal amplitude, but also its shape, 
possibly changing best-fit concentrations and slopes. 

The properties of these distributions, and the results 
of performing our fits to the averaged signals, are found in 
Table Q The important lesson is that the averaging pro- 
cess has little effect on the best-fit concentration parameter, 
as one would hope. The average signal can be modeled as 
an NFW profile to a high degree of accuracy. We remind 
the reader, however, that the process we have described in- 
corporates only the averaging over a distribution in halo 
masses, but not any scatter in the mass-luminosity rela- 
ti onship or the concentration -mass relationship. As shown 
in lMandelbaum et al ]|2005b'), the scatter in these relation- 
ships can have more profound efi'ects on the signal than the 
simple average over distributions. 

Another averaging process that may concern us is the 
average over redshifts. Because the nonlinear mass changes 
with redshift, the predicted value of concentration parame- 
ter c also changes with redshift. Hence, this may cause some 
blurring of the profile that will affect best-fit parameters. 
However, since our average over halo masses is also an aver- 
age over concentration parameters, and was found to have 
little effect on the profile shape, and because the change of 
concentration with z is small, this effect is also not expected 
to be significant. 

2.3 Method 

Our method is to fit physically-motivated models with the 
sum of a stellar plus dark matter component to AE. In 



practice, these 6-parameter fits, with four parameters for 
the dark matter and two for the baryonic profile, are highly 
degenerate due to the scarcity of information about small 
scales, so in order to obtain results, we make physically- 
motivated assumptions. For each sample, we do the follow- 
ing fits in order of increasing complexity: 

(i) Fix to a one-parameter fit, a singular isothermal 
sphere (SIS), and fit only for the mass. 

(ii) Fix to a more general single power-law profile, fitting 
for the mass and the logarithmic slope. 

(iii) Fix a, (3, and the stellar component; fit for c and 
NFW mass. 

(iv) Again, fix to an NFW profile, but try fitting for other 
profile parameters such as inner or outer logarithmic slope. 

One fact that may help us in these fits is that AE to 
some degree mixes information about p at different scales. 
For given transverse separations R, AE(_R) includes infor- 
mation from lower values of r due to the averaging E(< R). 



3 DATA 



The data used here are obtained from the SDSS (lYork. et alJ 
t200Ct). an ongoing survey to image roughly tt steradi- 
ans of the sky, and follow up approxima tely one million 
of the detected objects spectroscopically ijEise nstein et alJ 
l2nnil: lR.ichards et"alll2003l : IStraxiss et alJl2002l ). The imag- 
ing is carried out by drift-scanning the sky in ph oto- 
metric conditions LHogg c t al. 200 1; Ivc zic ct a l]| 2004l). in 
five bands (ugriz) ijFukugita et al.lll99fil: ISmith et al.ll2002h 
using a specially-designed wide-field cam^ra ijounnetejl 
1998). These imaging data are the source of the Large- 
Scale Structure (LSS) sample that we use in this pa- 
per. In addition, objects are targeted for spectroscopy us- 
ing these data jBlanton et alJ l2003j) and are observed 
with a 640 - fiber spectrograph on the same telescope 
jGxinn et alJ l200(t) . All of these data are processed by 
completely automated pipelines that detect and mea- 
sure photometric properti es of objects, and astrometri- 
cally calibrate the dat a (jLuoton et al.l 1200 ll : IPier et al] 
I2OO3I: iTucker et aljlioosl) . The SDSS is we ll underway, and 



has had six major d ata rel eases ( Stoughtgi^^^ 
Abazaiian ct al. 2003, 2004, 200^; iFinkbeiner et al 
Adelman-McCarthv et al..2006j . 



3.1 Lenses 

The galaxies used as lenses are those targeted as 
the spectroscopic Lum inous Red Galaxy (LRG) sample 
jEisenstein et al.l200J^ . including area beyond Data Release 
4 (DR4). The total area coverage is 5154 square degrees, as 
available in the most rece nt version of the NY U Value Added 
Galaxy Catalog rVAGG. iBlanton et aLllioOsI ). 

We include these galaxies in the redshift range 0.15 < 
z < 0.35, where the upper cutoff is designed to ensure 
that the lenses still have sufficient number of sources be- 
hind them. Within these redshift limits, the sample is 
approximately volume-limited with a number density of 

1.2 X 10~''(/i/Mpc)^, once we apply the additional color- 
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Table 1. Approximate values of minimum, maximum, mean, and median predicted masses for the two LRG samples using the halo mass 
function dn/dM for two cosmologies, and the best-fit masses to the model signal. 

Sample M^i„ M^ax (M) M^^^a-a (c) M„/ 

w,fit ^nfw.fit 

IQI^/i-IMq lO^^h-'^MQ W^^h-'^MQ lO^h-iM© IO^/i-IMq 
flm = 0.3, as = 0.9 

Faint LRGs 2.5 6.3 3.8 3.5 7.4 3.8 7.3 

Bright LRGs 6.3 oo 14 9.8 6.3 13 6.1 

n„ = 0.25, as = 0.75 

Faint LRGs 1.6 3.2 2.2 2.0 6.5 2.2 6.2 

Bright LRGs 3.2 oo 7.1 4.9 5.5 6.0 5.0 



magnitude cuts described in Loh & Strauss (2006, in prep.)^ 
to eliminate contamination from fainter, bluer galaxies be- 
low z ~ 0.23. (We note that this is the number density 
before we make our cut requiring that the galaxy be a host 
galaxy; the numbers in Table|5|give a lower average number 
density.) 

Model (not Petrosian) magnitudes were used for all lu- 
minosity cuts described in this paper; in all cases, h = 1 
magnitudes were used, so the division into samples does not 
depend on h. Magnitudes are corrected for extinction using 
reddening maps from SchleEel et al, L1998i) . We apply a k-l-e- 
correction (combined k-correction and correction for evolu- 
tion of the spectrum) to all magnitudes to redshift z ero as in 
Wake , et. al. (2006, in prep.) using the Bruzual fc Charlol 
i2003h stellar population synthesis code"^. We note that if 
we do not use these k-|-e-corrections , but simply use KCOR- 
RECT v4_l_4 (iBlanton et al.ll2003dl to the sample median 
redshift of 0.27, the resulting absolute magnitude histogram 
suggests that the galaxies above the median redshift are on 
average brighter than those below it. This result is consis- 
tent with the sign of passive evolution, for which galaxies are 
brighter in the past. In order to avoid our luminosity cuts 
selecting intrinsically different samples at high and low red- 
shift, we apply all luminosity cuts using the k-|-e-corrected 
luminosities. Figure |3 shows information about the redshift 
and magnitude distribution of the sample. 

In order to study variation with luminosity, we divide 
the sample into two subsamples at Mr = —22.3; roughly 2/3 
of the sample is fainter than this value, and 1/3 is brighter, 
giving approximately equal signal to noise in each bin. Sam- 
ple parameters are included in Table |5| 

Scale radii for the baryonic components were obtained 
by finding the average de Vaucouleurs profile fit radius for 
each sample from Photo, the SDSS processing pipeline; all 
radii quoted here are half-light radii. 

For this paper, we isolate "host" versus "satellite galax- 
ies" using spectroscopic galaxy counts in cylinders of comov- 
ing radius 2 /i~^Mpc and line-of-sight length Av — ±1200 
km s~^. Since we only want to eliminate galaxies for which 
there is another, brighter spectroscopic LRG nearby, we sim- 
ply require that those LRGs in our "host" sample either (a) 
be the only one in the cylinder, or (b) be the brightest in 
the cylinder. 

We thank Yeong-Shang Loh for providing the necessary files to 
implement this cut. 

^ We thank David Wake for providing the k+e-corrections for the 
appropriate models. 



LRG sample magnitude, redshift distributions 




Figure 3. Upper left panel: scatter plot of absolute magnitude in 
r band, k-corrected to z = 0.3, versus redshift; line shows average 
trend with redshift (total change of 0.2 magnitudes across the 
full redshift range). Upper right: absolute magnitude distribution 
dn/dMr when k-|-e-corrected to ^ = 0, with a line indicating 
the cut value between the faint and bright bins. Lower right: 
Mr k-|-e-corrected to z = versus redshift; line shows average 
trend with redshift, which is nearly constant (total change of 0.02 
magnitudes across the full redshift range). Lower left: comoving 
number density as a function of redshift, with a line indicating 
the average value. 



A problem with this scheme is fiber collisions; since 
fibers cannot be placed closer than 55", corresponding to 
200 h~^kpc at these redshifts, we may in principle still have 
some satellite galaxies that will distort the observed profile 
shape. Note that roughly 9 per cent of spectroscopic LRGs 
lack redshifts due to either fiber collisions, or due to redshift 
failure. We address this issue by using all galaxies passing 
the spectroscopic LRG sample cut I (below z = 0.36) and 
the additional stringent color-magnitude cut from Loh & 
Strauss (2006, in prep.), whether or not a spectrum was ob- 
tained, in order to test whether a galaxy is a host or a satel- 
lite. For those LRGs without spectra, we used the redshift 
of the nearest LRG if there was one within 55"; otherwise. 
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we used photometric redshifts f rom Kphotoz v4_1_ 4 with 
luminous red galaxy templates jBlanton et al ] |2003bl) . with 
az ~ 0.03 for these bright red galaxies. If the galaxy with- 
out a spectrum turned out to be brighter than a galaxy with 
a spectrum within a cylinder, we use the brighter one de- 
spite the lack of a spectrum. This cut on local environment 
reduces the sample by roughly 8 per cent. 

We obtained stellar mass e stimates in the following way . 
First, we use the estimates from lPadmanabhan et aP i2004) . 
which found that M stellar / L ~ 2>Mq/Lq for ellipticals (this 
is also con sistent with the relationship between Msteiiar and 
L seen in iMandelbaum et alJ l2005dl for bright red galax- 
ies). The average luminosities of the faint and bright LRG 
samples determined in terms of Lq using the absolute mag- 
nitude of stellar luminosity determined in iBlanton et alJ 
hOOSd) are 5.2 and 8.Qh~^ x lO^^ie respectively, and with 
L* = 1.2h~^ X W^'^Lq, we find (with appropriate factors of 
h = 0.7) stellar masses of 2.2 and 3.7 x lO"/i"^M0, respec- 
tively. 

Finally, we will further split the "bright" LRG sample 
at Mr — —22.6 in order to better trace the M{L) relation. 
When this split is done, the Mr < —22.6 sample contains 
1/9 of the fuU LRG sample, and the -22.3 > Mr > -22.6 
sample contains 2/9. 



3.2 Sources 

The source sample used is the same as that originally de- 
scribed in Mandclbaum et al. (2005a), hereinafter M05. This 
source sample includes over 30 million galaxies from the 
SDSS imaging data with r-band model magnitude brighter 
than 21.8, with shape measurements obtained using the 
REGLENS pipeline, including PSF correction done via re- 
Gaussianization (Hirata & Scljak 2003) and with cuts de- 
signed to avoid various shear calibration biases. A full de- 
scription of this pipeline can be found in M05; here we in- 
clude only a brief summary. 

The REGLENS pipeline obtains galaxy images in t he r 
and i filters from the SDSS "atlas images" (Stoughton c t alJ 
|2002^ . The basic principle of shear measurement using these 
images is to fit a Gaussian profile with elliptical isophotes 
to the image, and define the components of the ellipticity 



1 - (b/'^ f 
1 + {h/aY 



(cos 2(j>, sin 2<j)), 



(16) 



where b/a is the axis ratio and (j) is the position angle of 
the major axis. The ellipticity is then an estimator for the 
shear, 



(7+,7x) = 



(17) 



where TZ ~ 0.87 is called the "shear responsivity" and repre- 
sents the response of the ellipticity (Ea. ll6ll to a small shear 
CKaiser et al. 1995: Ber nstein & Jarvis 2002). In practice, a 
number of corrections need to be applied to obtain the el- 
lipticity. The most important of these is the correction for 
the sme aring and circularization of the galactic images by 
the PSF: IMandelbaum et aP Jiooi^) uses the PSF maps ob- 
taine d from stellar images by the psp pipeline iLupton et alJ 
120011) . and corrects for these using the re-Gaussianization 
technique of iHirata fc SeliaP JioO^, which includes correc- 
tions for non-Gaussianity of both the galaxy profile and the 



PSF. In order that these corrections can be successful, we re- 
quire that the galaxy be well-resolved compared to the PSF 
in both r and i bands (the only ones used for shape mea- 
surement), where we define the Gaussian resolution factor 



R2 = l 



(18) 



where the T values are the traces of the adaptive covariance 
matrices, and the superscripts indicate whether they are of 
the PSF or of the galaxy image. 

M05 includes a lengthy discussion of shear calibration 
biases in this catalog; we will only summarize these is- 
sues briefiy here. Our source sample is divided into three 
subsampl es: r < 2L r > 21, and high-redshift LRGs 
ySiscnstcin et al.ll200j) . defined using color and magnitude 
cuts as in M05 using s electi on criteria related to those from 
lEisenstein et al.H20o J) andlPadmanabhan et alJ^l2005^ . Us- 
ing simulations from lHirata^^enam2003l) to estimate the 
PSF dilution correction and analytical models for selection 
biases and other issues that affect shear calibration, we place 
2(7 limits on the shear calibration bias of [—0.05,0.12] for 
r < 21, [-0.08, 0.18] for r > 21, and [-0.06, 0.19] for LRGs. 

As shown in Eq. 13 the lensing signal AE is a product 
of the shear and factors involving lens and source redshifts. 
Since the lenses have spectroscopic redshifts, the primary 
difficulty is determining the source redshift distribution. We 
take three approaches, all described in detail in M05. For 
the r < 21 sources, we use photometric redshifts and their 
error distributions determined using a sample of galaxies 



in th e Groth strip with redshifts from DEEP2 (Davi^^^ 



20031: iMadgwick erai] 120031 ICoil et al.l Eooi 



Davis et al 



2005|), and require Zs > zi + 0.1 to avoid contamination 



from physically-associated lens-source pairs. For the r > 21 
sources, we use redshift distributions from DEEP2. For the 
high-redshift LRGs, we use photometric redshifts and their 
error distributions determined using data from the 2dF- 
Sloan LR G and Quasar Su rvey (2SLAQ), and presented in 
IPadmana bhan et al.l (|200^. Note that the LRG lenses used 
in this paper are typically at higher redshift than the lenses 
used in M05 and for this reason we are more sensitive to er- 
rors in the source redshift distribution. We will discuss the 
implications of this in H4.2I 

Finally, we have placed constraints on other issues af- 
fecting the calibration of the lensing signal, such as the 
sky subtraction problem, intrinsic alignments, magnification 
bias, star-galaxy separation, and seeing-dependent system- 
atics. As shown in M05 the calibration of the signal using 
the three source samples agrees to within 10 per cent, with 
a total la calibration uncertainty estimated at 7 per cent 
(r < 21) or 10 per cent (r > 21 and LRG). 



3.3 Signal computation 

Here we describe the computation of the lensing signal. Lens- 
source pairs are assigned weights according to the error on 
the shape measurement via 



Wis 



ai + a 



(19) 



SN 
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Table 2. The lens samples used for this analysis, including redshift, absolute magnitude, stellar mass limits, and typical scale radii. 



Sample 


N 


z limits 






A'/,, limits 


(L/Lq) 




^deV 














h-2l0l0 


(10"h-iMQ) 




Faint LRGs 


27 700 


[0.15,0.35] 


0.27 


0.24 


Mr > -22.3 


5.2 


2.2 


8 


Bright LRGs 


15 635 


[0.15,0.35] 


0.27 


0.24 


Mr ^ -22.3 


8.6 


3.7 


11 



where (t|^, the intrinsic shape noise, was determined as a 
function of magnitude in M05, figure 3. The factor of EJ^ 
downweights pairs that are close in redshift. 

Once we have computed these weights, we compute the 
lensing signal in 53 logarithmic radial bins from 20 /i~^kpc 
to 4 /i~^Mpc as a summation over lens-source pairs via: 



AE(i?) = 



(20) 



where the factor of 2 arises due to our definition of ellipticity. 

There are several additional procedures that must be 
done when computing the signal (for more detail, see M05). 
First, the signal computed around random points must be 
subtracted from the signal around real lenses to eliminate 
contributions from systematic shear. Second, the signal must 
be boosted, i.e. multiplied by B{R) = n{R)/nrand{R), the 
ratio of the number density of sources relative to the num- 
ber around random points, in order to account for dilution 
by sources that are physically associated with lenses, and 
therefore not lensed. 

In order to determine errors on the lensing signal, we 
divide the survey area into 200 bootstrap subregions, and 
generate 2500 bootstrap-resampled datasets. 

We also would like to account for the total calibration 
uncertainty of roughly 8 per cent at the la level (including 
redshift uncertainties, shear calibration, stellar contamina- 
tion, and all other effects that cause calibration bias). In 
principle, this can be done by multiplying the signal in each 
bootstrap resampled dataset by l-fCOSg where y is a Gaus- 
sian random number with mean zero and variance 1, and the 
net effect will be to cause correlations of order 8 per cent 
between radial bins in addition to any correlations due to 
shape noise, systematic shear, or other causes. However, we 
find that this procedure causes a significant bias in our fits, 
because points that are naturally below the best-fit signal 
have their errorbars increased by a smaller amount than 
those that are above the best-fit signal, and therefore those 
that are below the real signal are weighted more highly in 
the fits. Consequently, we generate the bootstrap-resampled 
datasets without including the calibration uncertainty, and 
add the 8 per cent uncertainty in quadrature with the sta- 
tistical error on the mass determination. 

Furthermore, to decrease noise in the covariance matri- 
ces due to the bootstrap, we rebin the signal into 18 radial 
bins total (of which 12 are in the range of radii that we will 
eventually elect to use for the fits). 

3.4 Theoretical systematics 

There are several theoretical systematics that are particu- 
larly problematic on small scales near massive galaxies. Here 
we describe two of them, estimate their magnitude, and ex- 
plain how we correct for them. 



0.01 




R [h-'Mpc] 



Figure 4. k, 74, and gt for the typical density profile described 
in the text, as labelled on the plot. 



3.4-1 Non-weak shear 

The first small scale systematic we consider is the fact that 
for these relatively massive halos, the assumption of weak 
shear may no longer hold on the smallest scales that we 
probe (20 /i~^kpc). In reality, what is measured in g-g lens- 
ing measurements is the reduced shear. 



9t = 



It 



(21) 



(where k = E/Ec). Hence, gt « 7t to a high degree of accu- 
racy in the 7t <c; 1, ft ^ 1 regime that is typically probed 
with weak lensing. Fig. 0| includes a plot of 7t, k, and gt for 
a typical LRG profile with the same density profile param- 
eters as in Fig.0 without adiabatic contraction (AC makes 
no perceptible difference in these quantities on these scales). 
As shown there, at the minimum scale at which we attempt 
to measure the lensing signal (20 /i~^kpc) the difference be- 
tween 7 and g is ~ 15 per cent, and at 40 /i~^kpc it is 9 per 
cent, causing us to overestimate AE. 

Unfortunately, it is difficult to take into account the 
effect of non-weak shear on our shear estimator, {e+), in 
a way that is completely model-independent. Appendix 1X1 
shows the rather involved calculations necessary to estimate 
the importance of any terms beyond first-order in 7t and 
K. As we show there (Eq. IA19I reproduced here for con- 
venience), the second-order calculation of the value of the 
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shear estimator is 

(e+) = 27^(AE)(E-l) (l + , (22) 

which gives a correction larger than the naive 1 + k. 

The two fractions in the correction term need to be 
estimated in some manner. The first one can be estimated 
using our earlier model of averaging AS over the halo mass 
function from the relevant minimum mass to the maximum 
mass, and is found to be « 1 for the faint LRG sample, 
and ~ 1.03 for the bright LRG sample. The second fraction 
in the correction term, (E~'^)/(E~^), can be very simply 
calculated when the signal is calculated, where the averages 
are computed using the weighting used to get the lensing 
signal, and typically equals approximately 1.05. Then, when 
the fitting is completed, we multiply the model AE by the 
correction term. 

3.4-2 Magnification bias 

The next issue we consider is the effect of magnification 
bias on the boost factors. Magnification bias will tend to 
increase or decrease the number density of observed sources 
relative to that expected from the random catalogs; hence, 
our procedure of multiplying the signal by n{R) /nrand{R) in 
order to correct for dilution by physically associated (non- 
lensed) sources will tend to overestimate the signal, since 
some of those additional sources are actually lensed. The 
degree of this overestimation was computed for each source 
sample as a function of its apparent magnitude and res- 
olution factor distributions, and was found in M05 to be 
5n/n = 1.9k, 0.7ft;, 3.1ft: for r < 21, r > 21, and LRG sources 
respectively. (Note that (E^^) varies with source sample as 
well, which gives different n for each sample.) In order to 
correct for this overestimation, we multiply the signal for 
each model by (1 + aE{E~^)) before comparing against the 
real data, where a and (E^^) are found as a weighted aver- 
age over the values for each source sample. Fig.^Jshows, for 
the model density profile used throughout this section, the 
value of It from which the size of this correction factor can 
be ascertained. 

A related issue is that of the effect of magnification 
bias on the source redshift distribution. This effect is not 
a problem for the r < 21 or the LRG source sample, since 
we use photometric redshifts that should still be reliable, 
so it is only a concern for the r > 21 sample for which we 
use a distribution p{z). Since magnification bias allows the 
observation of sources that would not have otherwise been 
seen, it means that in regions where magnification bias is 
significant, we tend to underestimate the source redshifts, 
and therefore E~^, overestimating the signal AE. The mean 
redshift of the r > 21 sample is (z) — 0.45. 

For the model LRG signal used in this section, we esti- 
mate the size of this effect as follows. We consider a mod- 
erately large value of «; = 0.1 such as may be present for 
transverse separations of 40 /i~^kpc (see Fig. |1J. For the 
r > 21 sample, this means that 5n/n 0.07, so the fraction 
of sources used that have been scattered into the sample 
is 0.07/1.07 ~ 0.065. We consider what happens if those 
sources are actually all at z — 1 (more than twice the mean 
redshift of this sample, giving a conservative estimate). EJ^ 
for these sources is then actually 50 per cent larger than 



what we have assumed. If 6.5 per cent of the sample has 
had its EJ^ value so underestimated, this means that the 
sample average E~^ value is thus 3 per cent too small, and 
the signal has been overestimated by this amount. Since this 
value is far smaller than other sources of systematic and sta- 
tistical error, and is a conservative estimate, we do not apply 
a correction. 



4 RESULTS 

Those who are uninterested in the details of the lensing sys- 
tematics tests in the following subsections may wish to skip 
directly to the results of the fitting, in subsection 14.41 

4.1 Small-scale software-related systematics 

One systematic potentially affecting all source samples on 
small scales is the sky subtraction syst ematic, an error in 
the sky estimation near bright lenses jMandelbaum et alJ 
l2005al : lAdelman-McCarthv et alJl2006^ . A similar system- 
atic is the possibility of errors in deblending. To test for 
these problems, we use simulations that were used for the 
small-scale LRG aut ocorrelation function measurement in 
iMasiedi et al.l i2005h . These simulations were created by 
placing fake LRGs on real SDSS images in run 2662, camcol 
1. For that paper, the issue of concern was the effect on a 
measured LRG flux due to the presence of another nearby 
LRG, so two fake LRGs were placed in each field at vari- 
ous angular separations. For this work, similar simulations 
were created with only one fake LRG per SDSS field, with 
470 fields total. The simulated LRG flux is r = 17.9, corre- 
sponding to a Mr = -22.5 LRG at 2 = 0.3 with r^ev = 3.4" . 
The axis ratio b/a — 0.7, with quantized random orientation 
relative to the scan direction (i.e., it was allowed to have any 
of 360 orientations evenly spaced in position angle). We used 
these values of flux and size to place conservative constraints 
on software- induced systematics; studying the variation with 
LRG size and flux is unfortunately prohibitive due to com- 
puter time and storage space limitations. 

The new images were then processed with the same ver- 
sion of the SDSS processing software. Photo, as the original 
data (rerun 137). We use these results to decide which ranges 
of transverse separations to use in the real data. The results 
presented in this paper are averaged over five separate "re- 
runs" with different simulated LRG positions but the same 
photometric properties. 

There are several concerns for which we would like to 

test: 

(i) Over- or under-estimation of the sky can lead to mod- 
ulation of the number density near LRGs (because they 
change the flux and apparent size, and hence scatter galax- 
ies into or out of our catalog). Since we assume that any 
change in number density relative to that around random 
points is due to physically-associated sources, this problem 
may cause a bias in the lensing signal. 

(ii) If the deblender or sky determination algorithm in- 
cludes a spurious gradient in the flux around LRGs, it can 
lead to the galaxies near the LRG (whether physically- 
associated or not) having some spurious additive and/or 
multiplicative shear. 
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(iii) If the LRGs cause nearby galaxies to have mis- 
measured fluxes, then depending on the way this effect varies 
with band, it may change the colors and hence lead to in- 
accurate photometric redshift determination. This problem 
would affect the r < 21 and LRG source samples. Further- 
more, for the r > 2f sample, for which we assume dn/dz, the 
redshifts of the sources may be underestimated, and hence 
AE overestimated, by some amount. 

In order to test for these effects, we searched for real 
galaxies around the positions of the simulated LRGs in the 
real data and the simulations. The first test, for modula- 
tion of number density, was performed simply by compar- 
ing the number of real galaxies around those locations in 
the simulation versus in the real data; results are shown in 
the upper left panel of Fig. |3 As shown there, there is the 
same depletion in number density on 30-90" scales as noted 
in previous works, with a crossover point occuring between 
20-30", and an excess of galaxies below that point. These 
findings are consistent with the sky being underestimated 
(and hence source galaxy fluxes overestimated) below the 
crossover point, and overestimated above that point. We 
note that the crossover point around 23" corresponds to 
~ 93 comoving /i~^kpc at the sample mean redshift. Also, 
in this panel as for the others, the results in the simulations 
match those in the real data above fOO". 

Next, we calculated both the tangential and 45 degree 
shear components around the positions of the fake LRGs, 
using only those real galaxies that passed cuts to be in our 
source catalog in both the real data and the simulation. 
The importance of this cut is that if we use the same real 
galaxies in both samples, then when we take the difference 
between the shears in these samples, all shape noise will be 
eliminated, giving an estimate of the "systematic" shear in 
both components. The resulting plot is shown in the lower 
right panel of Fig.|S] The x shear component, while contain- 
ing several points that are slightly discrepant from zero, is 
nonetheless consistent with zero spurious shear for all points 
above 10" ; below this scale, the errorbars on both shear com- 
ponents appear to be underestimated. The -I- shear compo- 
nent, the tangential shear, is also consistent with zero down 
to 10", or roughly 40 /i~^kpc: the constraint \'yaya\ < 0.004 
corresponds to |AEsys| < 20, which is far smaller than the 
statistical errors. Below this scale, constraints on contamina- 
tion are weaker, so we use a minimum transverse separation 
of 40 /i~^kpc for aU fits. 

This calculation only shows additive shear on galaxies 
that are matched. In the real data, any additive shear de- 
termined in this way will have to be multiplied by the boost 
factor B{R), since physically-associated sources are equally 
susceptible to software-induced systematics. Furthermore, 
this calculation only applies to galaxies that are matched 
in both the real and simulated data; it does not apply to 
galaxies that are scattered into the catalog due to sky un- 
derestimation. It is plausible that there is some selection 
bias operating in the selection as these galaxies, for instance 
causing us to select those that are aligned tangentially to 
the lens rather than radially, that would cause an additional 
bias in the results. Unfortunately, placing constraints on 
these galaxies is much more demanding because there are 
no matching galaxies in the real data with which they can 
be paired to eliminate shape noise. Finally, this calculation 



tells us about additive shear errors, but not about multi- 
plicative biases due to the deblender. 

The next question we address is that of the effect on the 
fluxes. While it is clear that either the fluxes or apparent 
sizes of sources near the fake LRGs is affected (since the 
number density of sources changes near them), the change in 
the fluxes may also affect classification into our three source 
samples, and may change photometric redshifts for the two 
samples that use them and the source redshift distribution 
for the r > 21 sample. For the galaxies that match in the 
real versus simulated data, the upper right panel of Fig. |K| 
shows the ratio of measured fiux in the simulated data versus 
in the real data, in both r and i bands. As shown, at 10" 
(« 40/i~^kpc), the fluxes are overestimated by roughly 15 
per cent, corresponding to a change in apparent magnitude 
of 0.15. Interestingly enough, the change is consistent for 
the r and i bands, even though the simulated LRG fluxes in 
these bands are not the same. This suggests that the effect 
may not noticably change photometric redshifts since they 
depend on color rather than magnitude. 

Finally, while we may guess that the change in number 
density is only due to the shift in the apparent magnitudes, 
the lower left panel of Fig. |K| indicates that the apparent 
sizes of the galaxies relative to the sizes of the PSF are also 
changed in a way that can cause the number densities to 
change as in the upper left panel. The quantity shown in 
this plot is 

sim real 

T(P) ^ T(P) ' ^ ' 

or the change in the trace of the adaptive moment matrix 
due to presence of the LRG, relativev to the trace of the 
adaptive moment matrix of the PSF. When we compare 
against Eq. 1181 we see that the change in resolution factor 
R2 due to the change in apparent size is 

5R2 = ^{1-R2)' (24) 

So, for a typical galaxy with R2 = 0.55 (sample median) 
at 10" from one of these LRGs, the resolution factor is in- 
creased by 0.09 to 0.64. A typical galaxy at separations of 
60" has resolution factor decreased by -0.01 to 0.54. 

Finally, we note for the sake of completeness that this 
simulation was only for one type of LRG, r — 17.9, with 
rdev = 3.3". This is reasonably close to the mean apparent 
magnitude of the spectroscopic LRGs used for this analysis. 
We have not derived variations of the effect with total LRG 
flux or apparent size, nor have we done any calculations for 
galaxies with profiles that are not purely de Vaucouleurs, so 
these results should not be used for, e.g., the Main galaxy 
sample which has a higher flux on average, and is a mixture 
of both de Vaucouleurs and exponential proflle galaxies. 

As a result of the flndings in this section that deblending 
and sky subtraction systematics become signiflcantly worse 
in a way that is difficult to quantify below about 10" (or 40 
/i~^kpc), in the analysis that follows, we only use the lens- 
ing signal above 40 /i~^kpc for fits even though plots show 
signal down to 20 h~^kpc . We correct the signal for the 
overestimation of the boosts due to the modulation of the 
number density near the LRGs, but because this correction 
is so large and uncertain below 10", we only show the re- 
sults below that separation to demonstrate consistency with 
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Figure 5. Four plots showing the effect of simulated LRGs on real galaxies in SDSS data. In all plots, the dotted line shows the nominal 
value if the simulated LRGs had no effects on the photometry of the real galaxies. The upper left panel shows the effect on the number 
density of real galaxies in our lensing source catalog around the positions of fake LRGs. The upper right panel shows the effect of the 
simulated LRG on the real galaxy fluxes in both r and i bands as labelled. The lower left panel shows the effect on the apparent size 
of the real galaxies relative to the size of the PSF, where the quantities measured are described in more detail in the text. Finally, the 
lower right panel shows the spurious shear caused by the deblender or the sky subtraction errors; in order to make the results easier to 
view, this quantity has been multiplied by (9/10")^. 



the results derived from higher angular separations. We also 
increase the error bars to account for systematics, with an 
increase of 50 per cent at 20 /i~^kpc down to no increase at 
50 ft-^kpc. 



4.2 Signal amplitude 

In this section we discuss the lensing signal calibration for 
each sour ce sample. The standard sy stematics tests as de- 
scribed in lMandelbaum et al.l ^2005a^ were performed to en- 
sure that systematics are under control (45-degree test, ran- 
dom catalog test), and no anomalies were found. 

One test that we focus on in particular here is the ratio 
test, which allows us to compare the effective average ampli- 
tude of AE for different source samples. The ratio t est is o f 
particular interest because, for iMand clbaum et alJ fcOOSaT) . 



the lensing signal for the three source samples was found to 
be consistent at the Icr level, but with lenses at a median 
redshift of 0.1. For this work, the lens redshift range is from 
0.15 to 0.35, with an effective mean of 0.24. The peak red- 
shift of the r < 21 sample is 0.35, quite close to the lenses, 
which means that for this work, we are far more sensitive 
to the d etails of the photometric re dshift error distribution 
than for ^ Mandelbaum et al.l ll2005all . In contrast, the mean 
redshift of the r > 21 sample is 0.45 and of the LRG source 
sample is 0.55, so we expect that the signal with these two 
source samples will still be fairly reliable. 

Fig- ISl shows the lensing signal for the bright LRG lens 
sample with the r < 21, r > 21, and LRG source sam- 
ples separately for 150 < R < 2000 /i~^kpc. We use these 
scales because they are the most free of systematics and offer 
the highest signal-to-noise. It seems clear that the LRG and 
r > 21 samples effectively agree in average amplitude, but 
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Lensing signal for bright LRG lenses 
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Figure 6. Lensing signal AS for the bright lens sample, with 
three different source samples as indicated on the plot. 



the r < 21 sample gives a lower overall calibration. When 
we determine an average effective amplitude of the lensing 
signal for each lens sample ("faint," "bright," and "all") for 
each bootstrap dataset, and compare this amplitude for the 
different source samples, we find the following relations hold: 



(AS),; 

(AE)l, 



= 0.98 ±0.11, 



(AS)r< 
(AE).; 



(AE) 



r>21 



(AE), 



(AE).>2i 
(AE)z,flG 



1.01 ±0.09, 



1.00 ±0.07, 



(AE) 



r<21 



(AE), 



(AE), 



(AE),; 



= 0.84 ±0.11 (faint) 
(25) 

= 0.86 ± 0.08 (bright) 
(26) 

= 0.85 ± 0.06 (all) 

(27) 

As shown in Eg. 1251 the signal for the LRG and r > 21 
source samples is statistically consistent for both LRG lens 
samples, but in contrast, the signal with the r < 21 (photoz) 
sources is about 15 per cent lower, a 2.5(t discrepancy when 
considering all lenses together. This finding suggests that for 
the r < 21 sample, some problem with the photometric red- 
shift error distributions is causing EJ^ to be overestimated 
by 15 per cent, thus underestimating AE by that amount. 

In order to test this hypothesis, we further split the 
r < 21 sources into two samples: those with photometric 
redshift lower than 0.4, and those with photometric redshift 
higher than 0.4. When we compare the signal with the "low 
redshift" and "high redshift" r < 21 sources, we find the 
following relations: 



(AE), 



,>0.4 



(AE), 



' 0.71 ± 0.13, Faint LRG lenses; 
0.80 ±0.12, Bright LRG lenses; (28) 
. 0.76 ± 0.09, All LRG lenses. 



The net result, then, is that the sources with photomet- 
ric redshift larger than 0.4 give signal that is 24 per cent 
lower than the signal with photometric redshift lower than 



0.4. When we take into account that this higher photo- 
metric redshift sample has more weight than the lower 
photometric redshift sample, this explains the signal with 
the full photometric redshift sample being low by 15 per 
cent. This result also explains o ur failure to notice t his dis - 
cre pancy in previous works; for [Mandelbaum et al] (jiooS^) 
and lMandelbaum et aP ll2005dl) . the mean lens redshift was 
lower, so (a) the sample with photometric redshift lower than 
0.4 had significantly higher weight than the sample with 
photometric redshift greater than 0.4; and (b) an imperfect 
understanding of the photometric redshift bias and scatter 
at photometric redshift larger than 0.4 gives a smaller frac- 
tional error on EJ^ for lower lens redshift. Consequently, 
the signal with the r < 21 sample was not as noticably dis- 
crepant in previous works. 

When considering E^^(zi = Zeff,Zs) as a function of 
source redshift, we note that for galaxies at r < 21, a net 
bias of 0.04 would be sufficient to achieve a 15 per cent 
change in overall calibration, e.g. if we assume that galaxies 
really at z = 0.4 are at z = 0.44. Alternatively, a smaller bias 
of ~ 0.02 and a modest increase in the net scatter of 0.02 
will tend to cause the same change in calibration. Since the 
peak of the redshift distribution for these galaxies is around 
z = 0.35, the photometric redshift error distribution was 
determined with a large number of galaxies from z — 0.3- 
0.4, but with a much smaller number of galaxies at z > 0.4, 
and hence a bias and scatter increase of these magnitudes is 
within the statistical error. 

Considering that we can fairly easily explain this bias 
in the r < 21 sample relative to the more reliable (for this 
purpose) r > 21 and LRG source samples, before doing the 
fits we have increased the signal with the r < 21 sample 
by 15 per cent as in Eq. 1251 and increased the statistical 
error by 6 per cent to account for the statistical error in 
the determination of this calibration factor. In practice, the 
increase in statistical error is achieved by multiplying the 
signal in individual bootstrap samples by a random number 
of mean 1 and Gaussian standard deviation 0.06. Note that 
this calibration change affects the mass determination, not 
the profile shape. 



4.3 Intrinsic alignments 

There are two intrinsic alignment effects that are impor- 
tant for lensing, but only one is important for galaxy-galaxy 
lensing as presented here. The intrinsic ellipticity-intrinsic 
ellipticity (II) correlation is not important for galaxy-galaxy 
lensing, because we average over all relative lens-source el- 
lipticity orientations. The second effect, which is important 
for this work, is an alignment between the intrinsic elliptic- 
ity and the local density field, or the GI correlation. This 
effect comes into play because we necessarily include some 
physically-associated pairs (i.e., pairs offenses and "sources" 
that are really part of the same local structure) , so if those 
sources have a tendency to align tangentially or radially rel- 
ative to the lens, they will provide an additive bias to the 
lensing signal. 

Correlations between intrinsic ellipticity and density 
have been dem onstrated robustly using bright, red galaxies 
on large scales jM andelbaum et al|2006l) and on small scales 
jMandelbaum et al.ll2005d : iDonoso et al.l 120061: lYang et al] 
l200^? ~with the effect being more controversial for the gen- 
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eral galaxy population and for spirals. In the cases consid- 
ered in these papers, the typical scenario is that a large, 
bright galaxy is found to point preferentially in the di- 
rection of nearby overdensities, e.g., along a cluster major 
axis. In our case, we are concerned that the fainter sources 
around a large, bright galaxy may be pointing preferen- 
tially radially or tangentially towards it; unfortunately, con- 
straints on this scenario are weak and /o r difficult to con- 
vert to measurements of l ensing shear jLee fc PenI l200lt 
iBernstein fc Norberd [20021: Iffirata et alJl2004^ . We expect 
that on scales larger than ~ 100 h~^kpc, above which the 
boost factors are fairly small, there will be little contribution 
from intrinsic alignments. 

Here, we try to measure the intrinsic alignment contam- 
ination of the lensing signal using photometric redshifts. For 
a given lens, the sources with r < 21 are split into a "phys- 
ically associated" (PA) sample with zi — 0.1 < Zs < zi + 0.1 
and a "lensed" (L) sample with Zs ^ zi +0.1. The difficulty, 
of course, is that the photometric redshifts have a broad er- 
ror distribution, so that both of these samples are going to 
have some galaxies in them that do not actually fit the cri- 
teria of being PA or L; in the two samples, the fractions of 
physically associated sources are f^f^^^ = (B*-^"*-' — 1)/B'^"*'' 
and f^p^2 = {B^^'> - l)/^'^'. We measure signals ■y*-^'^' and 
^f-^) = AE'-^^(E~^)<^' that are related to the "real" lensing 
and PA shears via 

^(P.) = + (1 _ /(f/))^f (29) 

- fp7W^/' + (1 - /r/^)AEi-)(E-)(--) 

where the lensing shear in the PA sample can be determined 
using AE from the lensing signal and multiplying by the av- 
erage E^^ for the PA sample. This average value can be 
determined by integrating over the known lens redshift dis- 
tribution, source photometric redshift distribution, and pho- 
tometric redshift error distribution, and is found to be a fac- 
tor of 2.5 lower than (E"^)'^', partly due to the significant 
number of sources in the PA sample that are foregrounds 
wth E^^ = 0. Thus, we can rearrange this equation to get 
T^M*'' intrinsic alignment shear in the PA sample, from 
the measured shear with the physically assocated sample. 

Next, we can estimate the contamination to the real 
lensing signal using the following equation; 



7 with PA sample 



0.01 



B 



(L) 



JpAnL + JpaIpa 



(30) 
(31) 



We note several things about this equation. The first term is 
simply the lensing signal, where the boost factor is necessary 
to account for the fact that not all the sources are lensed. 
The second term is the intrinsic alignment contamination 
term. We have pointed out with our notation a nuance that 
may be important: in Eg. I3UI we will be able to determine 
the intrinsic alignment shear due to physically associated 
sources in the PA sample; in Ea. l31l we actually need the in- 
trinsic alignment shear due to physically associated sources 
in the L sample. We return to this difference shortly. 

Fig. |7| shows the measured shear with the PA samples, 
.^(^^-4) ^jjg ^-^f^ jgjjg samples, and then the corrected 7p^'*'' 
after we remove the contamination due to lensing. As shown, 
the signals are consistent with zero above 100 /i~^kpc; below 




yC'l, bright LRGs 

ylP*), faint LRGs 

7<,P*>, bright LRGs 

.7™, faint LRGs 
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Figure 7. Measured shear with the PA source sample for both 
LRG lens samples, and the implied shear due to intrinsic align- 
ments after subtracting off the contribution due to lensing. 



that scale, they are consistently negative, indicating a ten- 
dency for satellites to "point" radially towards primaries. 
The for a fit to zero signal for the two lens samples 
are 11.6 and 12.5 for the bright and faint lens samples re- 
spectively (with 3 degrees of freedom with 40 < J? < 100 
/i~^kpc), giving a probability to exceed by chance of 0.01 
and 0.007, indicating a positive detection on these scales. If 
we average the value of Ip^"^' for those scales, we get aver- 
ages of -0.0084 ±0.0034 and -0.0086 ± 0.0030. We can also 
get the projected lensing contamination, if we make the cru- 
cial assumption that 7^^^' ~ by using (E^^) for the 
r < 21 sample of 1.9 x 10~*. This signal, when subtracted 
from the measured 7*^^^ to get a decontaminated lensing sig- 
nal, leads to an increase of the signal roughly equivalent to 
the size of the Icr errors or larger on these scales, and hence 
has a significant effect on fits to the profile. 

In comparison with recent works, we note that using 
SDSS data, jAEUstsson & Braincrd ( 20053 ^-l^o found a ten- 
dency for satellites to align radially towards primaries using 
spectroscopically-determined primary-satellite pairs. Over 
the range 7 < 7? < 50 /i~^kpc, they found a mean tangential 
shear of 74 = —0.030 ± 0.007. This number is, as expected, 
larger than our measured tangential shear, since ours is at 
a mean transverse separation that is close to three times as 
high as the mean separation for their measurement. 

However, we must consider further whether this is really 
the proper correction to apply. For the r < 21 sample, the 
sources that are truly physically associated in the L sample 
are, by definition, those with poor photometric redshifts. 
Hence, studies of the photometric re dshifts that we use 
in com parison with DEEP2 redshifts jMandelbaum et alJ 
l2005a) indicate that these will tend to be, on average, bluer 
and fainter. On the other hand, the physically associated 
sources in the PA sample are those with good photometric 
redshifts, which will tend to be brighter and redder. There 
is no a priori reason to suppose, then, that 7p^'*'' = 7pa, in 
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light of the fact that many intrinsic alignment measurements 
indicate that the effect depends sensitively on the galaxy 
luminosity and color. In fact, both large- and small-scale 
measurements seem to indicate that intrinsic alignments are 
negligible for bluer, fainter galaxies and large for bright, red 
ones, which suggests that 7^^] may well be zero, with no 
correction necessary l|Azzaro et a\ .l2006l:lMandelbaum et alJ 
l2006l:lYang et al.ll200gV 

The bottom line, then, is that it is not clear how to 
correct the signal with r < 21 sources for intrinsic alignment 
contamination in a way that is not based on many tenuous 
assumptions. We may come to a similar conclusion about 
the r > 21 sample. Since we only have a redshift distribution 
for this sample, the sample of galaxies that are physically 
associated are a mix of blue and red, but are on average 
(at a given redshift) fainter than the physically associated 
sources from the r < 21 sample, and hence the corrections 
derived from the r < 21 PA sample may not apply. 

Because of this problem with both the r < 21 and 
the r > 21 samples with respect to intrinsic alignments on 
R < 100 /!-~^kpc scales, we choose to rely solely on the high- 
redshift LRG sources on these scales. While the high-redshift 
LRG sample contains, by definition, bright red galaxies that 
may be expected to have strong intrinsic alignments, they 
are located at redshifts z > 0.4, whereas our lenses are in 
the range 0.15 < z < 0.35. Furthermore, we impose a cut 
on the source photometric redshift [zs > zi -\- 0.1) to fur- 
ther eliminate low-redshift contamination. The net effect is 
that the physically associated fraction at i? > 40 /i^^kpc is 
negligible when using this source sample, so it must be free 
of intrinsic alignments. To ensure that the signal on 40-100 
/i~^kpc scales, which may be crucial to understanding the 
profile, is as free of systematics as possible (while retaining 
higher statistical power on larger scales), we thus use only 
the LRG sources at i? < 100 ft~^kpc, and average all the 
source samples for R > 100 /i~^kpc. 

4.4 Fits to lensing signal 

We begin the fits with the simplest possible way, and add 
parameters until the best-fit does not decrease with the 
addition of parameters. Table|5]shows the results for best-fit 
parameters for the fits as enumerated and described there. 
In this table, in all cases the fits are from a minimum scale of 
40 /i~^kpc (to avoid systematics-dominated inner regions) to 
a maximum of 2 ^~^Mpc (to avoid the need for complicated 
modeling of the halo-halo term). Errors on all parameters 
are the formal Icr errors from the fit covariance matrix de- 
termined using derivatives of the with respect to each fit 
parameter. We remind the reader that the mass estimates 
have had an additional 8 per cent systematic error added 
in quadrature with the formal statistical error, and the fits 
are performed using the full covariance matrix (which in- 
cludes minimal correlations at small scales, and correlations 
as large as 20-30 per cent at large scales). 

4-4- 1 Faint LRG sample 

Accordingly, we begin by discussing the results for the faint 
LRG sample, with the simplest fit performed, that for a 
simple SIS with fixed stellar component, with only one free 
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Figure 8. Lensing signal for the faint LRG lens sample averaged 
over source sample, and best-fit signals. The vertical dotted line 
at 40 /i~^kpc shows the minimum scale used for the fits, and at 
800 /i~^kpc shows r^i^. The line types for each fit are indicated 
on the plot itself. We note that the signal shown here has been 
corrected for errors in the boost factor due to the sky subtraction 
error, but non-weak shear effects and magnification bias effecs on 
the boost factor are still incorporated in both the real and the 
model signal shown here. Errorbars are slightly correlated (20-30 
per cent level) on large scales. 



parameter. This fit (fit number 1 in Table |3J , along with 
several of the others to be discussed, appears in Fig. IHl with 
the measured signal. We have forced the profile in Eq. |S| 
to approximate a SIS by making very small, so that the 
logarithmic slope is equal to —jS — —2 for the entire radial 
range of interest. It appears from the top panel of Fig.|H|that 
the model overestimates the signal on 40-80 h~^kpc scales. 

Next, we consider the results for a general single power- 
law, fit 2. The addition of a single additional parameter, the 
exponent of the power- law (— /3), decreases the value by 
7.1 to 11.5 relative to fit 1, despite the fairly small change 
in exponent from —2 to —1.85 ± 0.06, suggesting that the 
addition of this fit parameter is justified. To properly evalu- 
ate the goodness of this fit relative to that for the SIS using 
the Ax^ value, we must take into account that it was calcu- 
lated using a b ootstrap cova r iance matrix, and as shown in 
Appendix D of iHirata et alJ (|20o3), the expected x^ values 
are not drawn from a x^ distribution in this case. The prob- 
ability value associated with exceeding this value of Ax^ by 
chance is p{> Ax^) = 0.04. On Fig. |S1 it is apparent that 
this model is preferred over the SIS because it does not over- 
estimate the signal on 40-80 /i~^kpc scales as much as the 
SIS (though the model signal is slightly high on 1-2 /i~^Mpc 
scales) . 

We then consider the simplest NFW-type fit, number 
3, which has the same number of fit parameters as the sin- 
gle power-law, and the same ■ Consequently, like the sin- 
gle power-law profile, the NFW model is preferred relative 
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Table 3. Fits performed, including identifying number, and fit parameters. Bold-faced fit parameters are meant to indicate that they 
were held fixed, whereas those in regular font are those that were allowed to vary, with la errors. For all fits, R^in = 0.04 ft-^Mpc. 
Note that the fitted values for concentration should be increased by about 20 per cent when compared to simulations to account for the 
scatter in cluster profiles, which was found to reduce the value in the lensing fits. Similarly, the fitted masses correspond better to the 
median rather than mean mass (Mandelbaum et al. 2005bl). 



Fit Description 


a 




c 


lO"/i-iM0 lO"/i-iM0 h-ikpc 


dof 





Faint LRG sample 



1 


SIS 


-2.0 


2.0 


10* 


2.4 ±0.3 


2.2 


8 


18.6 


10 


2 


Power-law 


-2.0 


1.85 ±0.06 


10* 


3.0 ±0.4 


2.2 


8 


11.5 


10 


3 


NFW 


-1.0 


3.0 


4.o±o.e 


i 2.9 ±0.4 


2.2 


8 


11.5 


10 


4 


NFW 


-0.50 ±0.18 


3.0 


7.7 


2.7 ±0.4 


2.2 


8 


13.8 


10 


5 


NFW 


1.0 


2.5 ±0.1 


7.7 


2.9 ±0.4 


2.2 


8 


10.1 


10 








Bright LRG 


sample 










6 


SIS 


-2.0 


2.0 


10* 


5.6 ±0.6 


3.7 


13 


28.6 


10 


7 


Power-law 


-2.0 


1.85 ±0.05 


10* 


7.8 ± 1.1 


3.7 


13 


16.7 


10 


8 


NFW 


-1.0 


3.0 


4.7±0.e 


; 6.7 ±0.8 


3.7 


13 


10.5 


10 


9 


NFW 


-0.72 ± 0.12 


3.0 


6.7 


6.4 ±0.8 


3.7 


13 


11.9 


10 


10 


NFW 


-1.0 


2.7±0.1 


6.7 


6.8 ±0.8 


3.7 


13 


9.4 


10 



to the SIS at the 96 per cent CL.* Visually, Fig. HI sug- 
gests that particularly on 40-80 /i~^kpc and 1-2 /i~^Mpc 
scales, the NFW model is preferred. The best-fit value 
of mass, (2.9 ± 0.4) x IO^^/i'^Mq. T his should be com- 
pared t o the median mass of the sample iMandelbaum et all 
l|2005a) . It is lower than the value expected naively from 
the halo mass function using high normalization, in Tabled 
it is close to the predicted lower mass limit of the sample, 
2.5 X li)^''^h~^ Mq, or to the expected mass value using the 
lower normalization model. (Actually placing constraints on 
Q.m and as using abundances is a more involved topic that 
we will address in a future publication.) 

We can also compare concentration parameter fits to 
theoretical expectations. Scatter in the mass-concentration 
relationship can lower the value in 2-d lensing fits when 
compared to average in 3-d. When using the signal for cen- 
tral galaxies in the brightest luminosit y bin in the simu- 
lations from IMandelbaum et alj (|20053), which incorporate 
both scatter in the mass-luminosity relationship and in the 
concentration-mass relationship, we find that the best-fit 
concentration is about 20 per cent lower than that expected 
for the corresponding mass. If we assume the same level of 
underestimate here, then this means that the fitted values 
should be increased by 20 per cent. We use the same value 
as an estimate of systematic error associated with this pro- 
cedure, but this source of systematic error could be reduced 
once large enough samples of simulated clusters are used to 
test the procedure. Thus the best-fit value of concentration 
is 5.0 ± 0.6(stat) ± l(sys), which should be compared to 7.7 
for the high normalization cosmological model and 6.5 for 
the low normalization model. We see that there is some ten- 



^ We have used the Ax^ values from a distribution assuming 
200 bootstrap regions, 12 radial bins, and 2 fit parameters. Since 
the SIS is not a special case of NFW, the actual distribution of 
Ax^ is more complicated in a way that is difficult to account for 
analytically, but it implies that the value of 96 per cent is overly 
conservative. 



sion with the high normalization cosmology, but the method 
needs to be tested further before we can draw conclusions 
about cosmology. 

For fit 4, we fix 13 to its canonical value of 3.0, and c 
to 7.7, but allow a and mass to vary. We find that a tends 
towards the maximum permitted value in the fits, —0.5, or a 
very shallow inner profile. This fit has a of 13.8, compared 
to fit 3, which had a fixed to the canonical value of —1.0 but 
allowed c to vary, with a of 11.5. This small difference 
in suggests that there is little difi'erence between these 
models, and that we are thus not highly sensitive to some 
combination of the values of c and a. 

For fit 5, we instead fix a to its canonical value of —1.0, 
and c to 7.7, but allow {3 and mass to vary. /3 tended to 
suggest a slightly shallower outer asymptotic slope, — 2.5 ± 
0.1 rather than —3.0, but still significantly different from the 
SIS value of —2.0. Again in comparison with fit 3, the 
for this fit is roughly 1 lower, which means this fit is slightly 
better, but not in a way that is statistically significant. 

Next, we tried several other fits, all of which had large 
degeneracies and thus none of which are shown in Table |21 
These fits include: 

(i) A fit for c, M, and Msteiiar' The best-fit value of 
Msteiiar was quite close to the assumed value, but with little 
constraining power. This is not too surprising, since our fits 
use scales Rmin ~ 5rdeV, well beyond the bulk of the stellar 
component. Nonetheless, it is useful to know that signifi- 
cantly different values of Msteiiar than the assumed value 
are not required by this fit. 

(ii) Simultaneous fits for NFW mass and either inner or 
outer slope along with concentration: we lack the power to 
constrain both the concentration and one of the slopes. As 
suggested by a comparison of fits 2-4 in Table |3 and Fig.|H| 
c and a are highly degenerate, with even fairly different 
combination of higher c and shallower a allowed with little 
change to the lensing signal. A comparison of fits 3 and 5 
shows that the same is true for /3 and c. 

Hence, while we are able to rule out the SIS model rel- 
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Figure 9. Lensing signal for tlie bright LRG lens sample aver- 
aged over source sample, and best-fit signals. The dotted line at 
40 /i~^kpc shows the minimum scale used for the fits, and at 
1 /i~^Mpc shows r^ir- The line types for each fit are indicated 
on the plot itself. Errorbars are slightly (20-30 per cent level) 
correlated on large scales. 

ative to NFW for this lens sample at the 96 per cent CL, 
we cannot place strong constraints on all profile parameters 
with this dataset. 

Bright LRG sample 

Next, we discuss the fits for the bright LRG sample, fits 6- 
10 in Table 01 with some fits shown with the real signal in 
Fig El The first fit we consider is the simple one-parameter 
fit (for M) with an assumed SIS profile and fixed stellar 
components, fit 6. The best-fit value of 28.6 gives a p- 
value of 0.004, suggesting that this model is a poor fit to 
the data From Fig. |51 it is apparent that the SIS model once 
again overpredicts the signal for 40-100 /i~^kpc scales. 

Fit 7 is for a single power-law profile, with the exponent 
allowed to vary from the SIS value. As for the faint LRGs, 
the best-fit exponent is — 1.85±0.05, with a large decrease in 
best-fit of 11.9 relative to fit 6, and a p(> Ax^) = 0.006, 
so the general power-law profile is a markedly improved fit 
over the SIS. Fig.j^lindicates that (as for the faint LRG sam- 
ple) this is the case primarily because this model predicts a 
lower signal than for the SIS model on 40-100 /i~^kpc scales, 
at the expense of overpredicting the signal on 1-2 /i~^Mpc 
scales. 

Fit 8 is for c and M with fixed NFW a, /?, and stel- 
lar components. The best-fit x^ is smaller than that of the 
SIS model (fit 6) by 18, and therefore the NFW model is 
favored. The SIS model in fit 6 is ruled out relative to the 
NFW model at the 99.96 per cent CL, and the power-law 
model (fit 7) relative to the NFW at the 94 per cent CL. 
As shown in Fig. |51 the NFW model does a better job pre- 
dicting the signal on 40-100 /i~^kpc and 1-2 /i~^Mpc scales 
than either of these other two models. This fit suggests a 



mass of 6.7 ±0.8 X lO"/i"^M0, so once again smaller than 
the predicted mean and median mass of this sample from 
the high normalization model halo mass function in Table Q 
and closer to the lower limit, or to the predicted mean and 
median halo mass using the lower o-g cosmology. 

For the concentration parameter in the NFW profile, 
we find c = 5.6±0.6(stat) ± l(sys) once we apply the 20 per 
cent upward correction. This should be compared to a pre- 
dicted c = 6.7 in the high normalization model and c = 5.5 
in low normalization model. While the high normalization 
model predicts somewhat too high values, although within 
the errors, the low normalization model agrees well with the 
observations. 

In fit 9, we constrain the concentration to the predicted 
value for this mass of 6.7, but allow a to vary. As shown, the 
best-fit value is somewhat shallower than the assumed value 
of —1.0, at roughly the 3cr level. However, the x^ value for 
this fit is worse than for fit 8 by roughly 1.4, so the lower 
value of concentration is needed to match the signal (but 
the difference between the two values of concentration is 
not statistically significant). 

In fit 10, we constrain the concentration to the predicted 
value of 6.7, but allow 13 to vary from canonical NFW value. 
As shown, the best-fit value of 2.7±0.1 is roughly 3cr from the 
canonical value, and the x^ value of 9.4 is actually slightly 
better than fit 18 which has (3 fixed to its canonical value 
and c free. This difference in x^ (of 1.1) is not actually sta- 
tistically significant, however. As for the faint LRG sample, 
we lack the power to simultaneously contrain concentration 
along with either the inner or outer slopes. 

For this sample, we also tried fitting simultaneously for 
c, NFW mass, and Msteiiar- As for the faint LRG sample, 
the best-fit mass was very close to the assumed value, but 
with very large errors, so we lack the power to constrain 
this parameter. Also as before, we note that the distinctions 
between the three NFW models shown on Fig.j^lis very small 
compared to the difference between the SIS, power- law, and 
NFW models. 

4.4-3 Mass to light relation 

In this section we also include the results of an attempt to 
determine the relationship between LRG r-band luminos- 
ity and mass with finer binning in luminosity. The bright 
bin has been split into two bins at Mr = —22.6. In this 
way, we are able to split this sample, which has an average 
mass corresponding to massive groups, into one sample with 
groups and the other with clusters. Results for the masses 
and concentrations as a function of luminosity appear in Ta- 
ble [l] and signal for the brightest bin is shown in Fig. 1101 
As shown by this figure, even considering the small size of 
this lens sample, we still have significant ability to constrain 
the profile and differentiate between NFW and SIS due to 
the high halo mass. When we split the sample further, how- 
ever, our ability to determine the profile to high precision 
degenerates, so we show no higher luminosity bins beyond 
this. 

Concentrations have, as in the previous section, been 
adjusted upwards by 20 per cent. Masses have not been ad- 
justed to account for scatter in the mass-luminosity rela- 
tionship, because the degree of adjustment that is necessary 
is not yet clear. In general, however, this adjustment would 
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Figure 10. Lensing signal for tlie Air < —22.6 lens sample aver- 
aged over source sample. 



Table 4. Masses and concentrations as a function of luminosity 
with finer binning in luminosity. 
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Figure 11. Measured NEW concentrations witli statistical and 
systematic errors added in quadratures, shown as a function of 
best-fit NFW mass. The lines shown on the plot are predictions 
for several different cosmologies. 
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the faint LRG sample, but this diflterence is not statisticaUy 
sig nificant, and can also be explained by an underestimation 



Mr > -22.3 
[-22.6, -22.3] 
Mr < -22.6 
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11.0 



2.9 ±0.4 
4.4 ±0.7 
12.9 ± 1.9 



5.0 ± 0.6(stat) ± 
6.8±l.l(stat) ± 
5.0±0.7(stat) ± 



have to lead to an increase of the best-fit masses, and hence 
it is clear that the mass is increasing quite rapidly with lu- 
minosity (roughly AI oc L^). Thus, the luminosity of the 
central galaxy can be used as a rough measure of the clus- 
ter mass, at least for masses below 2 X lO^Mo, although 
the (possibly significant) scatter in this relation cannot be 
established by our statistical approach. This approach to 
finding clusters is much simpler than other methods, such 
as the number of red galaxies (i.e. richness) or integrated 
light within a certain window in real and redshift space. 

Fig. llll shows the measured c as a function of halo mass, 
along with the predictions for several cosmological models; 
as shown, somewhat lower normalizations are preferred. We 
postpone the comparison between the measured and pre- 
dicted abundances to a future publication, since more work 
will be necessary to understand the effects of scatter in the 
M{L) relationship on the measured NFW masses. 



5 INTERPRETATION AND CONCLUSIONS 

While our fits only use scales above 40 h^^kpc due to con- 
cerns about lensing systematics, Figs.|H|and|^include scales 
down to 20 /i~^kpc the minimum scale probed. As shown 
in Fig.|Hl the signal below 40 /i~^kpc is fit somewhat better 
by the single power-law models than by the NFW profile for 



l(sys)of stellar mass rather than by any problem with the NFW 
l(sys)profile. Fig. |5| shows that the NFW model does a satisfac- 
l(sys)tory job describing the data even below 40 h~^kpc for the 
brighter lenses. For the faint lenses, several changes in pa- 
rameters could raise the signal on scales below 40 h~^kpc, 
such as an increase in stellar mass or a steeper inner slope for 
the DM profile. In both figures, we show the signal for the 
NFW fit without stellar component, and it is clear that the 
inclusion of the stellar component improves the agreement 
between model and observation below 50 h~^kpc scales. 

One might conclude, from Figs.|S]and|5] that it would be 
difficult to distinguish between inner DM profiles with any 
technique. However, there are several reasons why this is not 
the case. First, in using the lensing signal to distinguish be- 
tween these density profiles, we are necessarily susceptible 
to some parameter degeneracies that are not a problem for 
other methods, such as those against which we compare be- 
low. Second, other methods have a different dynamic range; 
in particular both strong lensing and kinematic studies can 
probe to far smaller scales (< 10 /i~^kpc). When combined 
with methods such as weak lensing that can constrain the 
profile on larger scales, in conjunction with constraints on 
the stellar profile, it seems conceivable that tight constraints 
could be placed on the inner dark matter profile. 

We note that the fact that the models are either compa- 
rable with or larger than the lensing signal on 1-2 h~^Mpc 
scales implies that we have successfully extracted a sample 
of LRGs that resides almost exclusively in host halos rather 
than as satellites sitting at outskirts of clusters for which 
the signal would be reduced. 

There are a number of conclusions that we can draw 
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from these fit results. First, the division into stellar and halo 
components allows us to fit the data very well. Adiabatic 
contraction does not seem to be necessary for our description 
on the scales of interest {R > 40 h~^kpc) because of the rel- 
atively low stellar to halo mass ratio, but the addition of the 
stellar component improves the fit on small scales over the 
simple NFW profile and there is some evidence that it may 
even be underestimated for the faint sample. We note that it 
is possible that AC occurred at a significantly earlier stage 
in the galaxy's formation, when M stellar / Mhaio was much 
larger. In that case, more significant AC on larger scales 
would have caused a steepening of the profile, which would 
then approximately be pres e rved through its later m ergers 
jBovlan-Kolchin et aLlboOSl iKazantzidis et alJl2006l) . This 
is another possible explanation for the excess over the NFW 
plus stellar mass model on 20 < i? < 40 /i^^kpc scales; how- 
ever, due to the uncertain systematics in this regime, we do 
not attempt to place constraints directly. Furthermore, we 
are unable to constrain the form of the stellar component; 
on the scales under consideration, a point mass works nearly 
as well as a Hernquist profile. 

We compare these results with those from other meth- 
ods that are used to determine the density profiles. Re- 
cent X-ray studies ha ve found that the profile is consistent 
with t he NFW mod el fVikhlinin ct al."'2005VFukaza wa et alJ 
I2OO6I : iHumphrev et al.. .200611 . though with indications for 
poor groups (consistent in mass with our "faint" sample) 
of so me excess over NFW on small scales dVikhlinin et al.l 
[200^, similar to our finding that the data is above the NFW 
model with stellar component in Fig. |H|for the fainter LRG 
sample. Thus, our results seem to be consistent with those 
from X-ray studies. 

Some kinematic and strong lensing studies (e.g., 
i Rusin &: Mai 1200 ll. iKoopmans fc Treull200a . iHuterer fc Mai 
I2OO4I. iKoopmans et aLlbOOfiTl have suggested a logarithmic 
slope ~ —2 on scales below 10 kpc, where optical obser- 
vations provide velocity dispersion information and strong 
lensing occurs. It is entirely possible that in that regime 
adiabatic contraction actually is necessary to understand 
the profile, and leads to the steepening of the profile com- 
pared to NFW. It is also possible that due to selection effects 
in the strong-lensing samples, they are not entirely compa- 
rable to our LRG-selected sample of groups and clusters. 
We note that we have central velocity dispersions for the 
LRG sample used in this analysis, which has been used to 
leajn about AC and dark matter inside the optical radius 
jPadmanabhan et al.ll20o3) . In the future it would be in- 
teresting to combine the weak lensing information with the 
kinematic information to learn more about the mass profile 
in the inner parts of the clusters. 

Studies of galaxy kinematics on larger scales, out to sev- 
eral Mpc, have found a strong preference for the NFW profile 
over a SIS p rofile o n the same scales as this study. These 
works fe.g.. Ivan dor Marcl et al.l I2OO0I: pi^iano fc Girardil 
l2003l:lRines et alJl200a : .Rines fc Diaferioll200d) tvpicallv fo- 
cus on several tens of clusters with a large number of spectro- 
scopic redshifts per cluster, and determine the profile using 
the cluster infall patterns. Generally the clusters for which 
this technique is used are in a higher mass range than those 
considered here, 10^* - li)^^h~^ Mq. It is useful to know that 
our results, which focus on the average profile of groups and 



clusters in a lower mass range, have been found to be true 
for individual higher mass clusters as well. 

We can see from the fits to the two samples sepa- 
rately that a characteristic halo mass for LRGs overall is 
~ 5 X lO^^/i~^M0. This corresponds to the transition regime 
between groups and low mass clusters. This is not too sur- 
prising, since the comoving number density of these objects 
(1.2 X 10~*(/i/Mpc)^) is an order of magnitude larger than 
for large clusters, and hence the typical LRG halo mass must 
be significantly smaller than that of massive clusters. If we 
adopt the definition that clusters have virial masses above 
IO^^'Mq then, given that the best fit mass of our bright 
sample corresponds to this value when we use h — 0.7, 
half of the halos in the bright sample can be defined as 
low mass clusters. This is also consistent with the findings 
in Loh & Strauss (2006, in prep.) that the typical LRG re- 
sides in groups, with some in the field, and some fraction 
in large clusters. We find that on average, brighter central 
LRGs reside in more massive halos, and the relation scales 
as M cc over the range we observe. 

When fitting to 2 ft" ^ Mpc, a profile that has a change 
in logarithmic slope with scale seems to provide the best fit 
to the data; the SIS profile is ruled out more robustly than 
a general power-law profile. Specifically, for the faint LRG 
sample, the SIS is ruled out relative to the NFW profile at 
the 96 per cent CL; for the bright LRG sample, the SIS is 
ruled out relative to NFW at the 99.96 per cent CL. 

Use of the the canonical slopes a = — 1 and P — 3 gives 
best-fit concentrations of 5.0 ± 0.6 ± 1 and 5.6 ± 0.5 ± 1 for 
faint and bright lens samples, respectively, where the first 
error is statistical and second systematic, since the scatter 
in concentration-mass relation leads to a bias in the fits and 
we do not have sufficiently large simulated samples of groups 
and clusters to model this bias accurately. By splitting the 
brighter bin further, we are able to obtain two different con- 
centration measurements for group and cluster-sized halos 
individually, as discussed in H4.4.3I The results are lower 
than the expected values from high normalization ACDM 
simulation, as shown in Fig. 1111 If we allow i}m = 0.25 and 
(78 = 0.75 in accordance with the recent WMAP 3-year re- 
sults (SpcrEcl ot al. 2006), the predicted values of c are low 
enough that the discrepancies are reduced to an acceptable 
level (below la when both statistical and systematic errors 
are included). Note that in principle, the statistical errors 
are sufficiently small that the method can provide constrain- 
ing power to differentiate among currently viable cosmologi- 
cal models. We refrain from making strong cosmological im- 
plications from this comparison because the method has not 
been calibrated yet in simulations over this halo mass range. 
When applied to lower masses below IO^^Mq it was found 
that the method underestima tes the true concentration pa- 
rameter by about 20 per cent llMandelbaum et"al ]l2005hD . a 
correction we have applied in the present analysis as well. 

In general, therefore, our results are in accordance with 
the predictions of ACDM N-body simulation predictions of 
halo profiles and statistically the measurements have very 
high signal to noise. More detailed analysis in the future 
will be able to differentiate among the cosmological models 
using both the measured cluster profiles and the abundance 
information. 
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APPENDIX A: EFFECTS OF NON-WEAK 
SHEAR ON SHEAR ESTIMATOR 

Here we compute the relation between the mean elliptic- 
ity and the disto rtion matrix. This relatio n has previously 
been examined bv lSchneider &: Seit j il995') for the case with 
Lu = 0; they derived an equation that is essentially equivalent 
to our Eq. 1A7II . although they did not provide the Taylor 
expansion in g. The analysis in this paper only makes use 
of the terms through second order in the weak lensing ap- 
proximation (i.e. second order in k. or in since this is 
adequate for the LRGs at tens of kpc radii where k ~ 0.1. 
We do however present the full formalism here, including 
all orders in in the hope that it will be useful in the 
future. 

The analysis here makes extensive use of the fact that 
our ellipticities transform covariantly under shear (this is in 
fac t true for all the adaptive ellipticities, e.g. as considered 
bv lBernstein fc Jarvid 12002?): note that many of the other 
ellip ticities used in the literature, such as "KSB" elliptici- 
ties jKaisereiaOlliii), do not share this property and the 
results of this section do not apply to them. 

The local relation between the source and image planes 
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(Al) 



is usually described by the Jacobian: 

J _ dxs _ f 1 — 7+ — ft — 7x+t<J 
dxi \ -7x - 1 + 7+ - K 

where the rows indicate the 1 and 2 components oi xs and 
the columns indicate those of a;/. Our basic problem here 
is to relate the mean ellipticity of the galaxies, (e), to the 
entries in J. In particular, we will find that (e) is not simply 
proportional to 7. 

It is most convenient here to work with the singular- 
value decomposition of J rather than with J itself. The 
singular-value decomposition is J = R((^)DR(S), where D 
is diagonal and R(q;) represents a counterclockwise rotation 
by angle a. If we define ^ = Tr D/2, then this can be written 
as 











R{e) = AR{ip)T{5)R{e). 



(A2) 

The four variables {7+,7x,k, oj} in the Jacobian have thus 
been replaced with the four variables {g,6,A,tp}. The ex- 
plicit formula for the conversion can be derived from the 
rectangular-to-polar conversion, Eq. (A5) of lHirata fc SeliakI 
i2003l) . which gives: 

2(1 -k) = 2Acos{ip + e), 
-2io = 2Asm{^p + 9), 

—27+ = —2gcos{Lp — 8), and 

-27x = -2gsm{^ - 8). 
The solution to this system of equations is 
A = ^{1-k)^+uj^, 



(A3) 



7+ + ll 

(1-k)2+Cj2 



— ( — arctan h arctan — ) , and 



1 / UJ 

— — arctan 

2\ 1 - K 



arctan 



7x 



(A4) 



The arguments of arctan in the last two lines have been 
written such that the numerator and denominator give the 
correct quadrant for the arctangent. Note that for the case 
of a single-screen lens with ui = 0, g reduces to the familiar 
"reduced shear" g = 7/(1 — k) (j^chncidcr & Scitz 1995). 

Of these four variables, A re-scales the source plane 
and so has no effect on the ellipticities. The rotation R(v') 
rotates the source plane and for randomly oriented source 
galaxies can have no effect on the average ellipticity. The 
transformation T{g) in Eq. (IA2II can produce net ellipticity 
in the -|- direction; and R{9) rotates the image plane by 
angle 9 and hence (e) by angle 29. Thus the observable (e) 
depends only on g and 9, and moreover the dependence on 
9 corresponds trivially to a rotation. 

Our problem is thus reduced to determining (e) for the 
special case where J — T{g). In this case symmetry consider- 
ations imply (fix ) = 0, so we need only determine (e+). If the 
intrinsic ellipticity of a galaxy is e'°\ then the equation for 
the fin al ellipticity is given by Eq. (2.5) of iMiralda-Escudel 
Jl99lfl: 



+ S 



5e. 



(0) ' 



(A5) 



w here the distortion 5 is rela ted to g by Eqs. (2-7) and (2-8) 
of iBernstein fc Jarvij (|2002|) : 



S = tanh(2 arctanh g) 



2<? 



1+5" 



(A6) 



It is simplest to evaluate (e+) for a population of sources 
whose intrinsic ellipticities have fixed magnitude e'"-* and 
random orientation <!>', other source populations can be han- 
dled by further averaging over the appropriate distribution 
for e'"-*. In this case, we may write e^' 



e'°' cos 2(f), so that 



1 

2^ 
1 

2^ 



2" cos2(j> + S 
l + 5e(0) cos 20 

27r 

-1 







+ 



d(f> 

5-^ 



1 + 5e(o) cos 2(l> 



d(j) 



5-r 



(A7) 



(The first term in the second line is integrated trivially, while 
the second term can be reduced to a rational integral by 
the substitution z — tan0.) We are interested in the weak- 
shear regime where we can use the Taylor expansion in 5, 
which is easily obtained from the binomial expansion of (1 — 

52g(0)2)-l/2. 



l_ie(0)2 
2 



+ 7:e 



(0)4 



16 



(0)2 



3g(0)4 

8 



+ Il6^ 



(0)6 



128 ' 



(A8) 



(There is no singularity at 5 = because the coefficients of 
cancel.) In the weak lensing literature, g is usually used 
instead of S; substitution of Eq. lA6t into Eq. HASII yields 



(e+) = ^ Cn 



2n + l 



(A9) 



where the first few coefficients in the expansion are 



Co 

Ci 

C2 

C3 



,(0)2 



2 - e 
-2 + 5e 
2 



(0)2 



-3e 



(0)4 



13e(")^ + 21e(»'* 



lOe 



(0)6 



and 



-2 + 25e' 



(0)2 



78e('"^ + 90e(°)« 



35e 



(o)( 



(AlO) 



For an ensemble of source galaxies with different intrinsic 
ellipticities, the values e'"-* " should be replaced by the ap- 
propriate moments (e^°' ") of the source galaxy distribution. 

In the case of a single-screen lens, it is convenient to 
expand {e+) in powers of the lensing strength EJ^, since 
both 7 and k are proportional to this quantity. In particular, 
we have 



7-^ 



: AE+ e: 



7x 



AEx E"\ and k = E EI 



(All) 

Here E is the surface density at the location of the image 
and AS is the surface density contrast defined in Fourier 
space by 



AE+(0±iAEx(0 =e 



zt2i arctan(Z2 / ^ 1 ) 



E(0. 



(A12) 



In the coordinate system where AS is in the -I- direction. 
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Eq. <A9ll then becomes 

This is turned into a power series in E^^ using the binomial 
theorem for negative powers, 

(l_,)-(.n+.)^^^(^2n+,j^,^ (A14) 

which gives: 



-3 



(A15) 

We now consohdate terms with the same exponent by defin- 
ing fc = 2n + j + f : 



fc=i 



L(fc-1)/2J 



n = 



xe;:*". (Af6) 

For the purposes of this work, we only use the fc = 1 
term in Ea. lA16l (or equivalently, the n = term in Ea. lA13ll . 
since higher order terms are 3rd-order or higher in n and 7. 
This lowest order term gives 

(e+) = (2-e(°)^)l^§^ (A17) 

we identify the (2 — e'"'^) factor with 2TZ, twice the shear 
responsivity. We then do the lowest order Taylor expansion 
of this equation, which tells us that 

(e+) = 27^{AEE7' + AEEE-2) (A18) 

When we do the usual approximation that the lenses are 
statistically identical, which is reasonable for these volume- 
limited samples, then the correction due to non-weak shear 
is 



(e,)^27^(AE)(E-)(l + ^^^g) 



(A19) 



Note that the factor in parenthesis is not the same as 
1 -I- K, for two reasons. First, (AEE)/{AE) > (E) since 
the massive objects get weighted more heavily. Second, 
(E-2)/{E-i) > (E-i). Both of these issues will tend to make 
reduced shear more important than the naive argument in 
the text that it depends on the size of 1 -I- k. 
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